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Abstract 

We obtain the initial phase space distribution after relativistic heavy ion collision by the CGC 
shattering method incorporating the uncertainty principle and solve the semi-classical Boltzmann 
equation which includes the gluon radiation processes. We present as a function of time the 
attenuation rate of high pr partons, which have transverse momenta over 6 GeV/c, in the medium 
which is formed after relativistic heavy ion collision. We calculate the elliptic flow as a function of 
an impact parameter, time and transverse momentum and also present the polar anisotropy, which 
gives the initial condition for color filamentation. 
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I. INTRODUCTION 

One of conclusions drown from the experiments at RHIC for the first period of experi- 
ment is that the strongly interacting Quark-Gluon Plasma is formed after Au-Au collision 
at maximum energy QjQ] ■ This conclusion, which is widely accepted by the community, is 
compelling based on several experimental results: The elliptic flow data, which is under- 
stood to be buildup at early stage of collision process, can be explained by hydrodynamical 
evolution of a plasma. And away-side jet quenching can be understood by the energy loss 
into the dense medium. 

It however is worth to reconsider the results in parton cascade simulation whether the 
data can be reproduced with no hydrodynamics. This is the purpose of our note. 

We describe briefly the Monte-Carlo simulation solving the semi-classical Boltzmann 
equations of motion in the next section. We obtain the initial phase-space distribution of 
partons in Section III. We presenr the results of our numerical solutions, namely, elliptic 
flow and polar anistropy in Section IV and the jet quenching in Section V. We summarize 
and conclude in Section VI. 



II. PARTON CASCADE CODE 



Although a system of partons evolves according to the field equation of quantum field 
theory, QCD, it is much more intuitive and rich in structure to use quantum transport 
equations^, 4] based on the field theory. It however is very difficult to solve the equations 
and interprete the functions which are defined on the spin and color index 



classical Boltzmann equation is used thus to capture t 



5 . The semi- 



re main effects [6|. To solve the 
partonic Boltzmann equation, the partonic cascade codes 3, have been developed. In the 
recent studies^, Q, the parton cascade code(PCC), which includes gg — > ggg channel in 
the secondary collision, has been developed and applied to the heavy ion collision. In the 
study, the minimum transverse momentum to produce minijets from the primary collision 
of two nuclei is set at 1 GeV/c and the small angle scatterings, sinO > pt/E where pt is 
the minimum momentum transfer set at 0.5 GeV/c, are included. Solving the Boltzmann 
transport equation with a given initial distribution by using a Monte Carlo simulation, the 
author presented the energy density, number density as well as momentum isotropicity of 
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partons. The detailed description to solve the equations can be found in 

We are obliged to mension a few flaw we have in solving the evolution equations: It 
has been known that the parton system can be very dense in early stage after relativistic 
heavy ion collision so that the partonic wave functions can overlap among them. This could 
invalidate the kinetic theory. On top of that we assume the particle is born just at the 
collision point with no formation time so that they can make a collsion just after collision. 
The processes we consider in our study are 

99 ^ 99, qq, (2.1 
gq gq, (2.2 

gq gq, (2.3 

q a q b <-> q c q d , (2.4 

qq <-> qq, (2.5 

q a q b q°q d , (2.6 

99 999- (2-7 

We however miss some of basic channels such as qg — *> q^, qq — > 77, qq — > (77, which could 
provide important information on the system. 

The cross sections for the processes up to the leading order (LO) are rewritten here for 



convenience 
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and 



(2.19) 



dt dt 
where s, i, u are Mandelstam variables. 

n 

We use the following approximation for the cross section of gg — > ggg process[l2j, 

da 99 ^ 999 9C A a 3 s q\ 



dq\dydk\ " 2 (ql + vl) 2 

Q(k±Xf — coshy)Q(y/s — k± coshy) 



(2.21) 



klV(kl + q{ + ^ D y-4kl q { 

We set the coupling constant to be 0.3 in our program. The minimum momentum transfer 
for the collision is 500 MeV/c and the total cross section for gg — > gg with this momentum 
cutoff, for example, is about 10/GeV 2 which is about 4 mb. We usually set the K-factor to 
2 or 4 to include the higher-order diagrams. 

Once the initial phase space is given we sequentially make a collision among the partons. 
To find the collision event relativistically, we only use the Coulomb-like part of QCD in- 
teraction with the assumption that the closest point between two partons is the collision 
event. 

III. INITIAL PHASE SPACE DISTRIBUTION 

The parton phase space distribution produced after relativistic heavy ion collision is one 
of most important indredients for the parton cascade simulation but also of most difficult 
subjets which is still in progress to its understanding. One of ways to obtain the distri- 
bution is to make elastic scatterings between the partons of projectile hadron or nuclei 



and those of target. T 
or minijets|l^. 



A target, inis, 



lis , so-called, factorization method, has been used to calculate jets 
There are many schemes to obtain the parton distribution of 
hadron or inu lcii l )! which is needed for the factorization method, for examples, we can 
fold the parton distribution of nucleon and the nucleon distribution of nucleus [la. 1 191]. This 
method only gives relatively high partons and soft partons can not be obtained so that it 
gives a reasonable distribution at LHC energy but does not give good one at RHIC energy. 
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b(fm) 


# of partons 





5400 


1 


5180 


3 


4150 


5 


2800 


7 


1570 


9 


630 



TABLE I: The number of parton produced after relativistic heavy ion collisions. 



A quite different way to obtain the initial parton information has been proposed by 
Krasnitz et al 
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211 ] based on CGC[22l.l23l]. namely, shattering of the color glass condensates 
of two colliding nuclei, which we will use in our study. 

The transverse momentum distribution of the produced partons after the shattering of 
two CGC is implicitly given by the eq. (24) of reference 20(. We can choose the transverse 
momentum using the Monte Carlo sampling. The rapidity distribution is in general flat at 
central rapidity region and fall off almost linearly to zero. We assume that the distribution 
is flat in rapidity between -2.5 and +2.5 and falls off to zero at ±5. The number of parton 
produced can be obtained by integration over transverse momentum and rapidity. This 
number however depends on the saturation momentum, Q s , which is introduced in CGC 
scheme and estimated to be lGeV < Q s < 2GeV at RHIC energy. We set Q s = 1.7GeV and 
the total number of initial partons is about 5400 for the head-on collision at RHIC energy. 
Since this gives only momentum distribution, i.e., transverse momentum and rapidity, we 
need to make a further assumption: The transversal position of produced parton can be 
choosen any point within the transverse overlap region of two colliding nuclei assuming the 
CGC is flat. On the other hand the longitudinal position will be chosen according to the 
uncentainty relation in between — 1/Ap z and +1/Ap z . And the time when the partons are 
born will be explained in detail later. 

Table 1 shows the number of partons produced after a relativistic heavy ion collision 
when the saturation momentum of color glass condensate is chosen to be 1.7 GeV/c. 

When the impact parameter is nonzero, we can calculate the number of partons and the 
transverse position by assuming the CGC plate is circular of the nucleus's radius so that 
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the total number of partons produced from non-central heavy ion collision is proportional to 
the overlap area relative to that of head-on collision. And we can choose the (transversal) 
position according to the overlap function, Eq. 3.1, where the probability density is given 
by 

P(f;b) = T AB (f;b) (3.1) 
= 4p£p^R\ - (r - b/2)^R% - (f + 6/2)2, 

Sampling the position according to this probability density can be performed using 
acceptance-rejectance method, namely we sample (x, y) within the allowed region randomly 
and calculate the probability density at that position to give P(x, y; b). We generate a ran- 
dom number r and compare to P(x, y; b) / P(0, 0; b) . If r is less than P(x, y; b) / P(0, 0; b) , we 
accept the (x,y) but if r is greater than that, we reject that and sample another positions. 
We next assume that those partons are born within the logitudinal uncertainty, Az > 1/Ap z , 
so that the z can be choosen between — Az and +Az. On the other hand, the initial parton's 
formation time will be given by 1/E even after the CGC shattering. The average energy of 
partons in rapidity \y\ < 1 is about 1.03 GeV and that in 1 < \y\ < 2 is about 2.1AGeV, 
and those of 2 < \y\ < 3 is 5.53GeV. Thus the partons of low energy has formation time 
0.1 — 0.2 fm/c while the high energy partons have much shorter formation time. 

Fig. [1] shows the rapidity distribution of the initial state sampled by Monte Carlo Method 
on head-on collision at RHIC energy. 

Since the distribution is not perfectly symmetric nor flat, we expect the initial elliptic flow, 
v 2, is not zero and fluctuate from set to set. This fluctuation shows up further analysis of jet 
quenching and elliptic flow. Note that we only sample 10 sets of initial distribution because 
of the limited CPU time. The total energy liberated from the CGC shattering is 34.5 TeV 
out of about 40 TeV of total available energy. 

Fig. [2] shows the transverse momentum distribution which can be compared with Eq.24 
ofref. [20] . 

The number of jets, which have the transverse momentum over 6 GeV in this study, for 
example, is about 45 for Au-Au head-on collision at RHIC energy when we use the shattering 
mechanism while the number reduces to about 10 for the factorization method. 

Fig. [3] shows the polar angle distribution of parton's momentum and shows the parton 
are strongly forward and backword oriented but relatively flat in central (transversal) 
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FIG. 1: Rapidity distribution of the test partons sampled by a Monte Carlo method, 
direction. 



IV. MOMENTUM ANISOTROPY 

The primary partons, which are produced directly from the colliding nuclei, have az- 
imuthal symmetry in momentum distribution since the collision cross section between par- 
tons has azimuthal independent. The partons will have numerous collisions among them- 
selves after they are born and could develop anisotropy in momentum distribution if there 
is spacial anisotropic. This azimuthal anisotropy can be extracted systematically by ex- 
panding the number density as a function of azimuthal angle with respect to the reaction 



plane [24, 



25] 



d 3 N 1 d 2 N 

E-- = — -=— [v + 2v x co8<f> + 2v 2 cos2(f) +...}, (4.1) 
dp- 3 7T dp^dy 

where v\ is the directed flow and v 2 the elliptic flow. When a parton momentum is known, 
the elliptic flow can be calculated by the equation, 

V2 = < 4^4 > (4-2) 

Px+Py 
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FIG. 2: pt distribution of the test partons sampled by a Monte Carlo method. 

where the bracket denotes the average over the partons and the impact parameter vector 
and collision axis defines the reaction plane. 

Fig. H] shows the elliptic flow as a function of transverse momentum for the parton system 
formed with Q s = 1.7 GeV, K = 4 and 6 = 5 fm. We include only those of partons with 
\y\ < 2. Since we average over 10 different sets of test partons, Fig. HJshows large fluctuation. 
We however can see that the overall number is close to those of elliptic flow per parton shown 
in Fig. 20 of referance 1]. 

Another issue raised recently is the polar anisotropy: One of reasons is that strong orien- 
tation of partons into the collision axis can cause the color filamentation and the momentum 



can be further isotropic because of the filamentation effects[26j. The other is that once the 
momenta of partons are isotropic, the hydrodynamic equation of motion can be applicable 
to describe the evolution even though the system is not (locally) thermalized 27j . This polar 
isotropy can be studied also with this cascade program. We first of all can define the polar 
anistropy as a function of ploar angle as follows: Consider the number density into a solid 
angle dfl, 
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FIG. 3: 6 distribution of the test parton's momenta sampled by a Monte Carlo method. 
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FIG. 4: Elliptic flow for Q s = 1.7 GeV and K = 4 at b = 5 fm. 



then we can define the normalized density such that 



n{0,<j>) = n{e,(t>)/N 7 



(4.4) 



where Nj< is the total number of partons. When the system is isotropic, n = 1/Air. We inte- 
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FIG. 5: Polar anisotropy for Q s = 1.7 GeV and K = 2 at b = fm and t = 3fm/c. 

grate over azimuthal angle, which is approximately isotropic, and can obtain the normalized 
anisotropy density e = n(rj) — 1/2 where r] = cosd. Fig[5] shows the polar anisotropy at t — 3 
fm/c. 

The figure shows that the system never have polar isotropy even in the central rapidity so 
that the condition for the filamentation can be met but that of applicability of hydrody- 
namics is not. 



V. JET QUENCHING 
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High energy parton can loose its energy by two mechanism 
first one is gluon radiation and the other is collisional energy loss. Study shows that the 
elastic scattering, which is major channel of partonic level, causes far less energy loss than 
the radiation. On the other hand, the system formed after relativistic heavy ion collision 
is rapidly expanding and those high energy partons of running to the collisional axis do 
not make any serious scattering but those of high transversal momentum have a chance 
to go through the expanding medium and lose substantial amount of energy by inelastic 



scatterings which could cause the gluon radiation. The parton cascade code(pcc), which we 
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FIG. 6: The number of jets as a function of time: formation time=0.1 fm, K=2, b=0 fm. 

are using, has the radiation channel, Ed. I2.211 and the available energy for this channel is 
high enough, several GeV, so that it is ideal tool to study the jet quenching. We simulate 
the given 10 sets of initial phase space distribution and average over the sets. The initial 
phase space from CGC shattering mechanism does not have back-to-back high energy 
parton. We thus calculate the number of high transvesal partons, pr > 6 GeV/c, as a 
function of time and obtain the attenuation rate. Fig. [6] shows the number of high pt 
partons with 0.1 fm/c formation time as a function of time. 

We can see that the number of jets increases in early stage of evolution, t < 1 fm/c, and 
decreases rapidly until 10 fm/c since then. This shows that those high energy partons, which 
trevel to longitudinal (collisional) direction, make collisions among them and produces many 
more high jets. These jets make further collisions with medium within 10 fm/c and lose 
energy. The total attenuation rate is estimated to about 1/3 or 1/4. 

Fig. [7] shows the pt distribution of jet. These high energy or high p^ partons eventually 
fragment and become relatively low energy or low p? hadrons. We hope to resolve the 
theoretical issue on the hadronization from a parton(s) and apply to this study. 
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FIG. 7: pt distribution of jets. 

VI. SUMMARY AND CONCLUSIONS 

We have shown that the elliptic flow can be reproduced from the kinetic theory with 
reasonable initial distribution and the polar anisotropy shows that the system is not 
isotropic enough to apply a hydrodynamics but the momenta of partons, especially high 
energy partons, are strongly oriented to the forward or backward direction which gives the 
color filamentation mechanism applicable. We have shown that the attenuation rate of high 
transverse momenta jets is about 3 or 4. But we cannot give any conclusive result from this 
since we expect the jet will fragmetate and hadronize substantially. 

We plan in near future to improve the statistics by using more sets of test partons and 
implement the hadronization mechanism which is recombination for low energy partons or 
fragmentation for high energy ones. It could bring us more detailed understanding on the 
relativistic heavy ion collisions and the fundamental signals of quark-gluon plasma. 
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